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We have developed a method for extracting the high-frequency noise spectral density of an rf- 
SQUID flux qubit from macroscopic resonant tunneling (MRT) rate measurements. The extracted 
noise spectral density is consistent with that of an ohmic environment up to frequencies ~ 4 GHz. 
We have also derived an expression for the MRT lineshape expected for a noise spectral density 
consisting of such a broadband ohmic component and an additional strongly peaked low-frequency 
component. This hybrid model provides an excellent fit to experimental data across a range of 
tunneling amplitudes and temperatures. 



Environmental noise is a critical concern in all ap- 
proaches to quantum computation, as it ultimately limits 
their feasibility. Understanding the origin of noise, the 
way it couples to qubits, and approaches to minimize 
or eliminate it will be key to any successful implemen- 
tation of a large-scale quantum computer. In this re- 
gard, experiments that use qubits as spectrometers to 
probe their environment are of particular relevance 
In this letter we describe an experimental approach that 
yields spectrosopic information concerning the environ- 
ment surrounding a superconducting flux qubit. The 
method described herein is generically applicable to any 
qubit whose potential energy landscape is bistable. 

Many experimental techniques have been employed to 
characterize noise in superconducting qubits. Rabi oscil- 
lation, spin-echo, and free induction decay experiments 
all provide measurements of decoherence times (2|-[^. 
These decoherence timescales are metrics of an aggre- 
gate response of the qubit dynamics to the degrees of 
freedom of the environment. The statistics of these envi- 
ronmental degrees of freedom can be characterized by a 
noise spectral density S{uj). Given a microscopic model 
of an environment, one can calculate S{uj) and, in turn, 
the decoherence timescales cited above. Inverting this 
process to infer S{uj) from decoherence times can yield 
model-dependent results that provide indirect probes of 
S{uj). On the other hand, at very low u one can obtain 
S{lj) through direct measurements of the slow variations 
of qubit parameters 0-0] • 

Quantum tunneling experiments with strong coupling 
to the environment represent an alternate means of quan- 
tifying S{ijj). In this case, S{lo) changes the observed 
rate of tunneling, thus providing a direct probe of 5 (a;). 
Such dissipative tunneling has been observed in a wide 
range of quantum systems, including superconducting 
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devices [lO|, nanomagnets single-electron tunnel 

junctions with resistive electrodes [12|, and carbon nan- 
otubes [13] • The dissipative environment of supercon- 
ducting flux qubits in particular is dominated by noise 
at low frequencies. Experiments measuring the rate of 



macroscopic resonant tunneling (MRT) of flux between 
the two lowest energy states have consistently yielded 
resonant tunneling peaks as a function of qubit bias that 
have a Gaussian lineshape near their maxima. Theoret- 
ical models assuming S{uj) is strongly peaked at small 
Lo naturally produce this Gaussian lineshape How- 



ever, experimental data show excess tunneling rates in 
the tails of the peaks that is not explained by such theo- 
retical models [4, [3] • In this work, we quantitatively 
show that the non-Gaussian tails in the MRT lineshapes 
can be attributed to components of S{uj) at high frequen- 
cies. We find that S{uj), as obtained for our rf-SQUID 
flux qubits, is well described by a broadband ohmic spec- 
trum plus a component that is strongly peaked at low lj. 

In a MRT experiment [i3,[l3|j one prepares a flux qubit 
in either the right or left well (|0), |1)) of the double- well 
flux potential of an rf SQUID. One then measures the 
rate of tunneling into the opposite well when the energy 
levels in the two wells are closely aligned (the process 
labelled To-j-i in Fig. [11(a)). For tunneling between the 
two lowest energy levels, one can map the system onto a 
two-state Hamiltonian of the form 



A(7^]/2-Q(7j2, 



(1) 



where ax,z are Pauli matrices, A is the tunneling ampli- 
tude, e is the energy bias between the wells, and Q is a 
noise operator that couples the qubit to its environment. 
The noise spectral density can be written as 



S{u) = / die^"*(Q(t)g(0)). 



(2) 



Hamiltonian ([T]) is valid as long as e, A<^ ?iWp, where Tiujp 
is the energy spacing to the next excited state within 
each well. In writing Hamiltonian ([T|), we have explicitly 
assumed that the environment couples an effective flux 
signal into the qubit body. 

In the limit of small tunneling amplitude, A ^ VF, 
where W is the noise magnitude defined in Eq. ([5]) be- 
low, coherent oscillations of flux between the wells are 
overdamped. One can then introduce a transition rate 
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FIG. 1: (a) Depiction of the double-well rf-SQUID po- 
tential and MRT processes discussed herein. Fq-h rep- 
resents interwell tunneling between the two lowest lying 
states. ro-»2 represents interwell tunneling between the low- 
est state in the initial well and the first excited state in 
the final well, (b) Schematic of the CCJJ rf-SQUID fiux 
qubit. External fiuxes applied to each closed loop indicated 
as G L,i?,ccjj,q,LT). 



from the left well to the right well, ro^i(e) = r(e), with 
the reverse rate ri_j.o(e) — r(— e). Under the general 
assumption that the noise possesses Gaussian statistics, 
the tunneling rate is given by 14, 18, 19]: 



r(,.|;/.e-/Sxp{/^ 



{hujf 



(3) 



If S{lo) has a dominant term that is strongly peaked at 
low frequencies, which we denote Slp{uj), Eq. ([3]) leads 
to a purely Gaussian lineshape offset by Cp from the qubit 
degeneracy point at e = [ij, |l7| : 



where 
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{V denotes principal value integration). For an envi- 
ronment in thermal equilibrium at temperature T and 
for noise strongly peaked at low frequencies such that 
w <g; ksT /Ti for all relevant lo in SLFii-^), t he p arameters 



W and ep are related via = W'^/2kBT []J,ll7|. Note 



that, since only cumulants of S'lf('-^) appear in Eq. (j4]) as 
given by Eq. ([5]) , MRT data cannot provide spectroscopic 
information regarding Slf{uj), or indeed any feature of 
S{uj) on frequency scales \uj\ < W/h. 

To account for the shape of an MRT peak away from 
its maximum, we allow for both a strongly peaked low- 
Lu component of the noise, 5'L_F(aj), and a broadband 
component of the noise, SnFii^)- Writing ^(a;) = 
SLF{i^)+SHF{i^) and expanding the exponential e~*'^* 



in Eq. ([3]) up to second order in w [1J|, we obtain the 



following expression for the flux tunneling rate: 
A2 
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Equation ([6]) implies that, in contrast to Slf{<jj), the 
shape of Shf{sjj) does affect the bias dependence of the 
tunneling rate r(e). Measurements of F(e) can there- 
fore be used to obtain spectroscopic information about 
Shf{^)- To do this, we perform two Fourier transforms. 
First, we define a function 



(7) 



Second, combining Eqs. © and ([7]), we express Shf{^) 
in terms of F{t): 
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liiF{t). 



(8) 



The environment of a flux qubit is typically in ther- 
mal equilibrium, as shown by previous MRT experiments 
17j . In this case, one can express the noise in terms of a 
spectral function J(w): 



SHF{co)^n'J{oj)/{l-e-^^'^), 



(9) 



where /? = l/ksT and J{uj) can be treated as an an- 
tisymmetric function of frequency, J{uj) = {Shf{oj) — 
SHF{-i^))/fi^- Jiuj) is related to F{t) through Eq. dH]): 



J{uj) = 2iuj'^ / dt sin Lut In F{t). 



(10) 



Using Eq. (fTO|). one can extract J(w) from measured 
r(e) versus e. Although J(a;) obtained in this way 
does not fully vanish at low frequencies, we assume that 
Slf{oj) ^ Shf{^^) for w < W/h, thus ensuring that 
the separation between the low- and high-frequency noise 
components underlying Eq. ([6]) is well defined. The 
MRT rate given by Eq. ^ is sensitive to Shf{^) in 
the noise correlator only in the time interval t < h/W. 
This means that the features of J{uj) on frequency scales 
LO < W/h cannot be resolved in F(e). In particular, J(w) 
for uj < W/h (where Slf{lo) ^ 0) is effectively obtained 
by smooth continuation from outside this range. 

An example functional form for J{uj) is that is fre- 
quently discussed is: 



J(w)=77a;|w/LJe|'e-l'^l/'^ 



(11) 



where r] is a dimensionless parameter characterizing the 
strength of the noise, s describes the noise frequency de- 
pendence, and ujc is a high-frequency cutoff [20|. We as- 
sume that hujc 3> ksT, W, e for our work. As we show be- 
low, J{uj) obtained from our flux qubits is consistent with 
an ohmic environment, for which s= 0. Using Eq. (jlip 
with s = and substituting into Eq. ([6]) yields 



r(e) = ^ / dre^(-^.)-W'^-V2 
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where = {it / I3ujc) '^^ A is the renormalized tunneling 
amplitude and Tc = {uJc)~^- For rj-^Air, the dependence 
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of Ar {~ A) on both T and uic is very weak. Moreover, 
the role of Tc in the integral is to remove the divergence 
and its exact value does not significantly affect the inte- 
gration result. The lineshape is therefore insensitive to 
Ldc for hujc » T, M^. Equations @ and are the two 
key predictions for the form of r(e) in the presence of 
S'l_f(w), with and without S'hf(w), respectively. 

Measurements were performed with a compound- 
compound Josephson junction (CCJJ) rf-SQUID flux 
qubit [lH . Figure [TJb) shows a schematic of this device. 
Static flux biases $2 and $^ are used to balance the crit- 
ical current of the left and right minor loops. A static flux 
bias $£rp allows one to adjust the inductance of the qubit. 
Time dependent flux biases $^ and ^ccjj permit control 
of the persistent current in the main loop 
energy A, and energy bias e= 2 

left and right wells of the potential, where $° is the qubit 
degeneracy point. Both A and |/^| are functions of $ccjj- 
Time dependent external biases and were 

delivered with high-precision room temperature current 
sources through flux bias lines that were inductively cou- 
pled to the qubit loop and the CCJJ loop, each with a 
mutual inductance ~ 2 pH. Cold filtering on all fiux bias 
lines limited the available bandwidth to 5 MHz. 

The circuit was manufactured on a silicon wafer with 
thermal oxide, Nb/Al/Al2 03/Nb trilayer junctions, and 
three additional Nb wiring layers insulated from one 
another with planarized, high-density plasma-enhanced 
chemical vapor deposited Si02. We mounted this chip 
inside an Al shield on the mixing chamber stage of a di- 
lution refrigerator with a minimum base temperature of 
21 mK. All qubit parameters were calibrated as described 
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elsewhere [15|. For our qubit, we extracted critical cur- 
rent /c= 3.38 ±0.01 fiA, inductance Lg=338±l pH, and 
capacitance C= 185 ± 5 fF. 

The basic MKT rate measurement technique has been 
described elsewhere [l3|- Figure [SJa) shows example rate 
measurements versus e. The lowest energy MRT peak is 
very well separated from the next lowest resonant tun- 
neling peak, which occurs at e/h ~ 13 GHz above the 
first peak (see process To^2 depicted in Fig. [T{a)). For 
the half-decade in r(e) near the peak, the Gaussian line- 
shape (jH) is a reasonable description of the data. To 
describe the resonant peak away from its maximum, we 
allowed for finite SHF{i^) and extracted J(a;) from the 
data using Eq. (fTO| . The results, plotted in Fig. EJb), 
agree very well with a straight line with slope 77=0.41 
up to |a;/27r|«4 GHz, the highest frequency for which 
we collected r(e). Linearity of J(cj) implies an ohmic 
environment. The theoretical lineshape calculated via 
Eq. ([T2|) using rj— 0.41 is in excellent agreement with the 
experimental data [solid curve in Fig. [Jfa)]. 

Given the success of an ohmic model for SnFi^^), we 
proceeded to directly fit a larger set of experimental data 
to Eq. p^ . Figure a) shows measurements of r(e), 
both Fq^i and Fi^Oi for three different values of ^^cjj 
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FIG. 2: (Color online) (a) Measured MRT rate as a function 
of e at T = 21 mK and "l>ccjj/'^'o = -0.6344. The dashed 
curve is calculated using Eq. @ with W/h = 0.47 GHz. The 
solid curve is calculated using Eq. (fT2l) with W/h — 0.47 GHz 
and 77 = 0.41. (b) J{uj) extracted from the experimental data 
using Eq. (fT0|) . The solid line is J{lj) = 0.41a;. 



($J^jj/$o = -0.6344, -0.6355 and -0.6365) with solid 
lines indicating fits to Eq. (IT2t . We obtained fit values of 
Ar/h = 3.15±0.07, 1.38±0.08 and 0.60±0.04 MHz from 
top to bottom, respectively. The extracted temperature 
was T — 21 ± 1 mK, in agreement with thermometry 
mounted on the mixing chamber of the dilution refriger- 
ator. The fit width W/h = 0.47 ± 0.02, 0.47 ± 0.02, and 
0.49 ± 0.02 GHz from top to bottom, respectively. We 
obtained = 0.41 ± 0.03,0.42 ± 0.03,0.425 ± 0.05 from 
top to bottom, respectively. 

We also performed MRT rate measurements at 
*^ccjj— ~ 0.6344 $0 for a range of chip temperatures from 
21 to 38 mK. Figure [31(b) shows example measurements 
and fits for three temperatures. We monitored the chip 
temperature via refrigerator thermometry and confirmed 
it by measuring the qubit transition width as described 
in jl51] . The temperatures extracted from the fit parame- 
ters used in Fig ^h) match those reported by thermom- 
etry and those obtained via qubit transition width mea- 
surements. Fit values of 77 and W were relatively insen- 
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FIG. 3: (Color online) Example measurements of MRT rate 
versus e for (a) three different barrier heights ($|f;,jj/$o = 



—0.6344, —0.6355, —0.6365 from top to bottom, respectively) 
at T = 21 mK, and (b) a single barrier height, ^^^^j/^o = 
—0.6344, and three temperatures. The hollow (solid) symbols 
are ro->i(ri_>o). The solid lines are a fit to Eq. (fT2]). 



sitive to chip temperature over the range probed in our 
experiments. We note that the T-independence of W im- 
plies a 1 /T-dependence of the linear response of the low 
frequency environment (typical, for example, of a pop- 
ulation of paramagnetic spins) [17]. In contrast, the T- 
independence of r] implies T-independence of the linear 
response of the high frequency environment, consistent 
with an ohmic environment. We observed a gradual in- 
crease of Ar/h with T: Ar/h = 3.15 ± 0.07, 3.45 ± 0.08, 
and 3.6 ± 0.1 MHz for 21 mK, 30 mK and 38 mK, re- 
spectively. Only a fraction of this increase in Ar/h with 
T can be accounted for by the renormalization of A in 
Eq. (|T2t . given the best fit value of rj. The rest may be 
due to the influence of the higher energy levels in the flux 
potential wells. 

The dimensionless parameter rj characterizes the am- 
plitude of the high frequency noise spectral density at a 
given ^'ccjjj therefore one particular persistent cur- 
rent I /Pi, for our devices. To move closer to a physical 
picture of the source of high frequency noise, we can re- 
late S'//f(w) to an effective flux noise S^{uj) by scaling 
the former by the persistent current at which the mea- 



surement was performed, Sis,{uj) = Shf{<^) / ■ 

The source of ohmic flux noise can be parameterized as 
an effective resistance Rg shunting the qubit junctions. 
The noise spectral density would then be 
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Using the measured qubit properties and the extracted 
rj w 0.42, we calculate an effective shunt resistance 



Rs 



-^4^ — i - 20 \Vt. 
Tirj 



(14) 



Another potential source of ohmic noise could be external 
qubit flux bias leads, each of which can be modeled as an 
impedance coupled to the qubit body with M = 2 pH. 
Considering the bias coupled to the qubit body and using 
measured qubit properties and rj w 0.42, we calculate 
an effective impedance 



Re(Zo) = 



1 rj. 



(15) 



Both of the ohmic sources hypothesized above pre- 
dict impedances that are at least an order of magnitude 
smaller than expected (independent junction measure- 
ments suggest Rs > 500 kft; we estimate Zq ^ 25 ft), 
making the ultimate source of the high frequency envi- 
ronment uncertain. Generally, the amplitude of the high 
frequency noise should depend strongly on the details of 
the qubit wiring, junction size, and strength of coupling 
to bias leads depending on its source. Future measure- 
ments of this amplitude for a variety of qubits with a 
range of wiring and junction sizes will allow us to probe 
the ultimate source of this noise. 

To summarize, we have developed and experimentally 
tested a method for extracting the high frequency noise 
spectral density S'hf('^) from MRT rate measurements 
on flux qubits. Our experimental data are consistent with 
an ohmic spectral density up to aj/27r = 4 GHz. We have 
derived a theoretical expression for the MRT lineshape 
that includes both low and high frequency noise compo- 
nents. The resulting model fits the experimental data 
very well. In particular, this model explains tunneling 
rate measurements away the resonant peak, where the 
model without high frequency noise fails. Our method 
allows further exploration of high frequency noise in de- 
vices via its dependence on qubit geometry and fabri- 
cation details. A systematic study of a range of qubit 
designs will aid in ultimately understanding the origin of 
high frequency noise in superconducting qubits. 

We acknowledge fruitful discussions with F. Cioata, P. 
Spear, E. Chappie, P. Chavez, C. Enderud, J. Hilton, C. 
Rich, G. Rose, M. Thom, S. Uchaikin, and B. Wilson. 



5 



[1] R. J. Schoclkopf ct al., ArXiv Condensed Matter e-prints 

(2002) , arXiv:coiid-mat /02 10247. 
[2] R. W. Simmonds et al., Phys. Rev. Lett. 93, 077003 

(2004). 

[3] J. M. Martinis et al., Phys. Rev. Lett. 95, 210503 (2005). 
[4] D. Bennett et al., Quantum Information Processing 8, 
217 (2009). 

[5] R. McDermott, IEEE Transactions on Applied Supercon- 
ductivity 19, 2 (2009). 
[6] F. Yoshihara et al., Phys. Rev. Lett. 97, 167001 (2006). 
[7] R. C. Bialczak et al., Phys. Rev. Lett. 99, 187006 (2007). 
[8] T. Lanting et al., Phys. Rev. B 79, 060509 (2009). 
[9] F. Yoshihara et al., Phys. Rev. B 81, 132502 (2010). 
[10] R. Rouse, S. Han, and J. E. Lukens, Phys. Rev. Lett. 75, 



1614 (1995). 

[11] L. Thomas et al.. Nature (London) 383, 145 (1996). 
[12] W. Zheng et al., Sohd State Communications 108, 839 
(1998). 

[13] Y. Bomze et al., Phys. Rev. B 79, 241402 (2009), 
1010.1527. 

[14] M. H. S. Amin and D. V. Averin, Phys. Rev. Lett. 100, 

197001 (2008). 
[15] R. Harris et al., Phys. Rev. B 81, 134510 (2010). 
[16] T. Lanting et a.. Phys. Rev. B 82, 060512 (2010). 
[17] R. Harris et al., Phys. Rev. Lett. 101, 117003 (2008). 
[18] M. H. S. Amin and F. Brito, Phys. Rev. B 80, 214302 

(2009). 

[19] D. V. Averin, J. R. Friedman, and J. E. Lukens, Phys. 

Rev. B 62, 11802 (2000). 
[20] A. J. Leggett et al.. Rev. Mod. Phys. 59, 1 (1987). 



